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Abstract 

We discuss the motivations and perspectives for the studies of the mesons of the (bc) 
family at LHCb. The description of production and decays at LHC energies is given in 
details. The event yields, detection efficiencies, and background conditions for several 
B c decay modes at LHCb are estimated. 
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1 Introduction 


The B c meson is the ground state of be system which in many respects is an intermediate 
between charmonium and bottonium systems. However because the B c mesons carry flavour, 
they provide a window for studying the heavy-quark dynamics that is very different from those 
provided by cc- and 66-quarkonia. 

The first observation of approximately 20 B c events in the B c —> J^lv decay mode by the 
CDF collaboration [1] demonstrates the possibility of the experimental study of the B c meson. 

The be states have rich spectroscopy of the orbital and angular-momentum exitations. 
Below the threshold of the decay into B — D pair, one can expect 16 extremely narrow states 
which cascadcly decay into the ground pseudoscalar state with mass of about 6.3 GeV by 
radiating photons and pion pairs. The annihilation decays can occure due to weak interactions 
only and, hence, are suppressed for excited levels. 

The production mechanism for be system differs in an essential way from that for bb system, 
because two heavy quark-antiquark pairs must be created in a collision. While the bb pair 
can be created in the parton processes qq, gg —> bb at the order of a|, the lowest order 
mechanism for the creating of bc system is at least of qq, gg —> (6c)6c, and gluon-gluon 
contribution dominates at Tevatron and LHC energies. At LHC with the luminosity of about 
C = 10 34 cm -2 s -1 one could expect around 5 x 10 10 B c events per year. At Tevatron energy 
the expected yield should be at least one order smaller. 

The weak decays of B c mesons are attractive due to presence of both channels: a) the 
6-quark decay with the c quark as a spectator, and b) the c-constituent decay with 6 as a 
spectator. In addition the weak annihilation contribution into decay channels is quite visible 
(around 10%). 

The dominant contribution into B c life-time ( tb c ~ 0.5 ps) comes from the c-quark decays 
(~ 70%) while the 6-quark decays and weak annihilation add about 20% and 10% respectively. 

The accurate measurement of the B c life-time can provide us with the information on both 
the masses of charm and beauty quarks and the normalization point of non-lcptonic weak 
Lagrangian in the B c decays. 

The experimental study of the semileptonic decays and the extraction of the decay form- 
factors can test the spin symmetry derived in NRQCD and HQET approaches. The measure¬ 
ment of the branching fractions for semileptonic and hadronic modes can provide an informa¬ 
tion about the parameters of weak Lagrangian and hadronic matrix elements determined by 
the non-perturbative effects due to quark confinement. 


2 The mass spectrum of the (bc) family 

The most accurate estimates of (6c) masses [2, 3] can be obtained in the framework of non- 
relativistic potential models based on the NRQCD expansion over both 1/rriQ and v re i —> 0 

[4]- .’ 

The uncertainty of evaluation is about 30 MeV. The reason is the following. The potential 
models [5] were justified for the well measured masses of charmonium and bottomonium. 
So, the potentials with various global behaviour, i.e. with the different r —» oo and r —» 0 
asymptotics, have the same form in the range of mean distances between the quarks in the 
heavy quarkonia at 0.2 < r < 1 fin [6]. The observed regularity in the distances between the 
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excitation levels are approximately flavor-independent. The latter is exact for the logarithmic 
potential (the Feynman-Hell-Mann theorem), where the average kinetic energy of quarks T 
is a constant value independent of the excitation numbers (the virial theorem) [7]. A slow 
dependence of the level distances on the reduced mass can be taken into account by the use 
of the Martin potential (power law: V(r) = A(r/ro) a + C, a<Cl) [8], wherein the predictions 
are in agreement with the QCD-motivated Buchmiiller-Tye potential with the account for the 
two-loop evolution of the coupling constant at short distances [9]. 

So, one gets the picture of (be) levels which is very close to the texture of charmonium and 
bottomonium. The difference is the jj-binding instead of the LS one. 

The spin-dependent perturbation of the potential includes the contribution of the effective 
one-gluon exchange (the vector part) as well as the scalar confining term [10]. 
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The model-dependent value of effective a s [3] can be extracted from the data on the 
splitting in the charmonium 


M(J/T) - M( Vc ) = a s -^\R(0)\ 2 « 117 MeV. 

9 m* 

We take into account the renormalization-group dependence of a s at the one-loop accuracy 
by means of introduction of the quarkonium scale [2] 

H 2 = (p 2 ) = 2 (T)m red . 


The estimated difference between the masses of basic pseudoscalar state and its vector exci¬ 
tation [2] is equal to 

M(B* + ) - M(B+) = 65 ± 15 MeV. 

The mass of the ground state [2] equals 

M(B+) = 6.25 ± 0.03 GeV, (2) 

which is in agreement with the CDF measurements M(B C ) = 6.4 ± 0.19 GeV [1]. 


2.1 Radiative transitions 

The bright feature of the (be) family is that there are no annihilation decay modes due to the 
strong interaction. So, the excitations, in a cascade way, decay into the ground state with the 
emission of photons and pion-pion pairs. 

The formulae for the El-transitions are slightly modified. 

T (nPj -> rtSi + 7) = ^ a em Qls t 2 (nP; nS) wj(nP ) , 


3 



7.5 



4P 




4D 

3S 





BD threshold 


3P 




3D 

2S 

r 



0 



2P 

? 



i + 



o + 

IS 

l' 



0 



state 

Martin 

BT 

1% 

6.253 

6.264 

l 1 ^ 

6.317 

6.337 

2% 

6.867 

6.856 

2 1 S 1 

6.902 

6.899 

2 3 P 0 

6.683 

6.700 

2 P 1+ 

6.717 

6.730 

2 P 1 ,+ 

6.729 

6.736 

2 3 P 2 

6.743 

6.747 

3'Po 

7.088 

7.108 

3 P 1+ 

7.113 

7.135 

3 P 1 ,+ 

7.124 

7.142 

3 3 P 2 

7.134 

7.153 

3P> 2" 

7.001 

7.009 

3 5 P 3 

7.007 

7.005 

3 3 Ph 

7.008 

7.012 

3P> 2'“ 

7.016 

7.012 


Figure 1: The mass spectrum of (be) with account for the spin-dependent splittings. 
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In eq.(3) 

(4) 


where Q c ^ are the electric charges of the quarks. For the B c meson with the parameters from 
the Martin potential, one gets Q e g = 0.41. wj(nL) is the probability that the spin S — 1 in 
the nL state. Sjj / are the statistical factors. The I(nL\nL') value is expressed through the 
radial wave functions, 

I(nL] nL') = \ J R nL (r)R nU (r)r 3 dr\ . (5) 
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For the dipole magnetic Ml-transitions one has 


16 

rfn'Si n'Sf + 7 ) = Y l4a w 3 (2/ + 1) A 2 , , 


(6) 


where 

A if = I Rns(r)R nS (r)j 0 (ujr/2)r 2 dr , 


and 
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Note, that in contrast to the i/j and T particles, the total width of the B* meson is equal to 
the width of its radiative decay into the R c ( 0 _ ) state. 

Thus, below the threshold of decay into the BD-pair the theory predicts the existence of 
16 narrow (be) states, which do not annihilate due to the strong interactions, but they have 
the cascade radiative transitions into the ground long-lived pseudoscalar state, the B+ meson. 


state 

Ttot, KeV 

dominant decay 
mode 

BR, % 

l^i 

0.06 

l% + 7 

100 

2% 

67.8 

l^O + 7T7T 

74 

2 1 S 1 

86.3 

RSi + 7T7T 

58 

2 1 Pq 

65.3 

l^i+7 
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2 P 1+ 
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l^i+7 

87 

2 P l'+ 

139.2 

l 1 5'o + 7 

94 

2 3 P 2 
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1^i+7 
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3'P 0 

44.8 

2%+7 

57 

3 P 1+ 

65.3 

2%+ 7 

49 

3 P 1 /+ 

92.8 

2 %+ 7 

63 

3 3 P 2 

71.6 

2%+7 

69 

3 D 2~ 

95.0 

2P 1+ + 7 

47 

3 5 D 3 

107.9 

2 3 P 2 + 7 

71 

3 s D 1 

155.4 

2^0 + 7 

51 

3D 2 

122.0 

2 P l'+ + 7 

38 


Table 1: The total widths of (6c)-states with Martin potential 


3 B c lifetime and inclusive decay rates 

The B c -meson decay processes can be subdivided into three classes: 

1) the 6-quark decay with the spectator c-quark, 

2) the c-quark decay with the spectator 6-quark and 

3) the annihilation channel B+ —■> l + ui(cs,us ), where l = e, jx, r. 
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In the 6 —> ccs decays one separates also the Pauli interference with the c-quark from the 
initial state. In accordance with the given classification, the total width is the sum over the 
partial widths 


T(B C -»■ X) = r(6 -»■ X) ± r(c -»• X) ± r(ann.) + r(PI). 

For the annihilation channel the r(ann.) width can be reliably estimated in the framework 
of inclusive approach, where one takes the sum of the leptonic and quark decay modes with 
account for the hard gluon corrections to the effective four-quark interaction of weak currents. 
These corrections result in the factor of Gq = 1.22 ± 0.04. The width is expressed through the 
leptonic constant of fs c ~ 400 MeV. This estimate of the quark-contribution does not depend 
on a hadronization model, since a large energy release of the order of the meson mass takes 
place. From the following expression, one can see that the contribution by light leptons and 
quarks can be neglected, 

r(ann.) = ]T ^\V bc \ 2 - m%/m 2 Bc ) 2 • C { , 

l=T,C 

where C T — 1 for the r + /z T -channel and C c = 3|V r cs | 2 o 2 for the cs-channel. 

As for the non-annihilation decays, in the approach of the Operator Product Expansion 
for the quark currents of weak decays [21], one takes into account the oq-corrections to the 
free quark decays and uses the quark-hadron duality for the final states. Then one considers 
the matrix element for the transition operator over the bound meson state. The latter allows 
one also to take into account the effects caused by the motion and virtuality of decaying quark 
inside the meson because of the interaction with the spectator. In this way the b —> ccs decay 
mode turns out to be suppressed almost completely due to the Pauli interference with the 
charm quark from the initial state. Besides, the c-quark decays with the spectator 6-quark 
are essentially suppressed in comparison with the free quark decays because of a large bound 
energy in the initial state. 


B c decay mode 

OPE, % 

PM, % 

SR, % 

6 —> cl + vi 

3.9 ± 1.0 

3.7 ±0.9 

2.9 ±0.3 

b — > cud 

16.2 ±4.1 

16.7 ±4.2 

13.1 ± 1.3 

Eb~>c 

25.0 ±6.2 

25.0 ±6.2 

19.6 ± 1.9 

c —> sl + Vi 

8.5 ±2.1 

10.1 ±2.5 

9.0 ±0.9 

c —> sud 

47.3 ± 11.8 

45.4 ± 11.4 

54.0 ±5.4 

Ec -> s 

64.3 ± 16.1 

65.6 ± 16.4 

72.0 ±7.2 

B+ —> t + v t 

2.9 ±0.7 

2.0 ±0.5 

1.8 ±0.2 

B+ -»• cs 

7.2 ± 1.8 

7.2 ± 1.8 

6.6 ±0.7 


Table 2: The branching ratios of the B c decay modes calculated in the framework of inclusive 
OPE approach, by summing up the exclusive modes in the potential model [13, 14] and 
according to the semi-inclusive estimates in the sum rules of QCD and NRQCD [24, 25]. 

In the framework of exclusive approach, it is necessary to sum up widths of different decay 
modes calculated in the potential models. While considering the semileptonic decays due to 
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the b —> cl + ui and c —> sl + vi transitions, one finds that the hadronic final states are practically 
saturated by the lightest bound lS'-state in the (cc)-system, i.e. by the r] c and J/tp particles, 
and the lS'-states in the (hs)-system, i.e. B s and B*, which can only enter the accessible 
energetic gap. 

Further, the b —> cud channel, for example, can be calculated through the given decay width 
of b —> cl + vi with account for the color factor and hard gluon corrections to the four-quark 
interaction. It can be also obtained as a sum over the widths of decays with the (ud)-system 
bound states. 

The results of calculation for the total B c width in the inclusive OPE and exclusive PM 
approaches give the values consistent with each other, if one takes into account the most 
significant uncertainty related to the choice of quark masses (especially for the charm quark), 
so that finally, we have 

t[B+] 0PEiPM = 0.55 ± 0.15 ps, (8) 

which agrees with the measured value of B c lifetime. 

The OPE estimates of inclusive decay rates agree with recent semi-inclusive calculations 
in the sum rules of QCD and NR.QCD [24, 25], where one assumed the saturation of hadronic 
final states by the ground levels in the cc and bs systems as well as the factorization allowing 
one to relate the scmileptonic and hadronic decay modes. The coulomb-like corrections in 
the heavy quarkonia states play an essential role in the B c decays and allow one to remove 
the disagreement between the estimates in sum rules and OPE. In contrast to OPE, where 
the basic uncertainty is given by the variation of heavy quark masses, these parameters are 
fixed by the two-point sum rules for bottomonia and charmonia, so that the accuracy of SR 
calculations for the total width of B c is determined by the choice of scale /i for the hadronic 
weak lagrangian in decays of charmed quark. We show this dependence in Fig. 2, where 
^ < n < m c and the dark shaded region corresponds to the scales preferred by data on the 
charmed meson lifetimes. 

Supposing the preferable choice of scale in the c —> s decays of B c to be equal to /. l 2 Bc ~ 
(0.85 GeV) 2 , putting ai(fj,B c ) = 1-20 and neglecting the contributions caused by nonzero a 2 
in the charmed quark decays [25], in the framework of semi-inclusive sum-rule calculations we 
predict 

t[B c ] sr = 0.48 ± 0.05 ps, (9) 

which agrees with the direct summation of exclusive channels calculated in the next sections. 
In Fig. 2 we show the exclusive estimate of lifetime, too. 

4 Semileptonic and leptonic modes 

4.1 Semileptonic decays 

The semileptonic decay rates are underestimated in the QCD SR. approach of ref. [22], because 
large coulomb-like corrections were not taken into account. The recent analysis of SR. in 
[23, 24, 25] decreased the uncertainty, so that the estimates agree with the calculations in the 
potential models. 

The absolute values of semileptonic widths are presented in Table 3 in comparison with 
the estimates obtained in potential models. 
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Figure 2: The B c lifetime calculated in QCD sum rules versus the scale of hadronic weak 
lagrangian in the decays of charmed quark. The wide shaded region shows the uncertainty 
of semi-inclusive estimates, the dark shaded region is the preferable choice as given by the 
lifetimes of charmed mesons. The dots represent the values in OPE approach taken from ref. 
[21]. The narrow shaded region represents the result obtained by summing up the exclusive 
channels with the variation of hadronic scale in the decays of beauty anti-quark in the range of 
1 < /ib < 5 GeV. The arrow points to the preferable prescription of /./ = 0.85 GeV as discussed 
in [25], 



In practice, the most constructive information is given by the J / V I ; mode, since this char- 
monium is clearly detected in experiments due to the pure lcptonic decays [1]. In addition 
to the investigation of various form factors and their dependence on the transfer squared, 
we would like to stress that the measurement of decay to the excited state of charmonium, 
i.e. -0', could answer the question on the reliability of QCD predictions for the decays to the 
excited states. We see that to the moment the finite energy sum rules predict the width of 
B+ —> ili'l + u decays in a reasonable agreement with the potential models. 

4.2 Leptonic decays 

The dominant leptonic decay of B c is given by the tu t mode (see Table 2). However, it has 
a low experimental efficiency of detection because of hadronic background in the r decays or 
a missing energy. Recently, in refs. [26] the enhancement of muon and electron channels in 
the radiative modes was studied. The additional photon allows one to remove the hclicity 
suppression for the leptonic decay of pseudoscalar particle, which leads, say, to the double 
increase of muonic mode. 


5 Non-leptonic modes 

In comparison with the inclusive non-leptonic widths, which can be estimated in the framework 
of quark-hadron duality (see Table 2), the calculations of exclusive modes usually involves 
the approximation of factorization [11], which, as expected, can be quite accurate for the 
B c , since the quark-gluon sea is suppressed in the heavy quarkonium. Thus, the important 
parameters are the factors Gq and a 2 in the non-leptonic weak lagrangian, which depend on 
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Mode 

r [is] 

r [is] 

r [17] 

r [is] 

r [19] 

r [20] 

Bi 

—> v c e + iy 

11 

11.1 

14.2 

14 

10.4 

8.6 

Bt 

—> r]cT + L> 

3.3 



3.8 


2.9 

Bt 

-»■ Vc e+l/ 

0.60 


0.73 


0.74 


Bt 

-»■ v' c t + v 

0.050 






Bt 

— > J/t>e + v 

28 

30.2 

34.4 

33 

16.5 

18 

Bt 

— > J/V , t + z/ 

7.0 



8.4 


5.0 

Bt 

— > ip'e + u 

1.94 


1.45 


3.1 


Bt 

— » 

0.17 






Bt 

— » D°e + iy 

0.059 

0.049 

0.094 

0.26 

0.026 


Bt 

-> D°t + u 

0.032 



0.14 



Bt 

-»• D*°e + u 

0.27 

0.192 

0.269 

0.49 

0.053 


Bt 

-> D*°t + u 

0.12 



0.27 



Bt 

— > B Q s e + v 

59 

14.3 

26.6 

29 

13.8 

15 

Bt 

-»• Bfe+v 

65 

50.4 

44.0 

37 

16.9 

34 

Bt 

— » B°e + u 

4.9 

1.14 

2.30 

2.1 



Bt 

-»• B*°e + u 

8.5 

3.53 

3.32 

2.3 




Table 3: Exclusive widths of semileptonic B+ decays, T in 10 15 GeV. 


the normalization point suitable for the B c decays. 

The QCD SR estimates for the non-leptonic decays of charmed quark in B c give the 
agreement of results with the values predicted by the potential models is rather good for the 
direct transitions with no permutation of colour lines, i.e. the class I processes with the factor 
of aq in the non-leptonic amplitude determined by the effective lagrangian. In contrast, the 
sum rule predictions are significantly enhanced in comparison with the values calculated in the 
potential models for the transitions with the colour permutation, i.e. for the class II processes 
with the factor of ci 2 - 

Further, for the transitions, wherein the Pauli interference is significantly involved, the 
class III processes, we find that the absolute values of different terms given by the squares 
of ai and aq calculated in the sum rules are in agreement with the estimates of potential 
models. However, we stress that we have found that clue to the Pauli interference determining 
the negative sign of two amplitudes with (i\ and <22 the overall sign in some modes should be 
different from those obtained in the potential models. Taking into account the negative value 
of a -2 with respect to aq, we see that half of decays should be enhanced in comparison with the 
case of Pauli interference switched off, while the other half is suppressed. The characteristic 
values of effects caused by the Pauli interference is presented in Table 4, where we put the 
widths in the form 


T = r 0 + AT, T 0 = x\ a\ + X 2 al, AT = za\ 02 - 

Then, we conclude that the Pauli interference can be straightforwardly tested in the listed 
decays, wherein its significance reaches about 50%. 

At large recoils as in B+ — > J/ x Tn + ( y p + ), the spectator picture of transition can be broken 
by the hard gluon exchanges [27]. The spin effects in such decays were studied in [28]. However, 
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Mode 

Ar/r 0 , % 

Bi - V c Df 

59 

Bt VcD* + 

-41 

Bt - J/^Df 

-55 

Bt -> j/^h; + 

53 

B+ - ricD + 

43 

Bt -»• VcD* + 

-47 

Bt - J/^D+ 

-46 

Bt -> J/^D*+ 

48 


Table 4: The effect of Pauli interference in the exclusive non-leptonic decay widths of the B c 
meson with the c-quark as spectator at a\ = 1.18, <22 = —0.22. 


we emphasize that the significant rates of B c decays to the P- and D-wave charmonium states 
point ont that the corrections in the second order of the heavy-quark velocity in the heavy 
quarkonia under study could be quite essential and suppress the corresponding decay rates, 
since the relative momentum of heavy quarks inside the quarkonium if different from zero 
should enhance the virtuality of gluon exchange, which suppresses the decay amplitudes. 

The widths of non-leptonic c-quark decays in the framework of the sum rule are greater 
than those of potential models. In this respect we check that our calculations are consistent 
with the inclusive ones. So, we sum up the calculated exclusive widths and estimate the total 
width of B c meson as shown in Fig. 2, which points to a good agreement of our calculations 
with those of OPE and semi-inclusive estimates. 

Another interesting point is the possibility to extract the factorization parameters oq and 
a 2 in the c-quark decays by measuring the branching ratios 


rft + - 

-> B+K° } 

T[B+ - 

-> B + K *°] 

U-B+ - 

-»• B°K+ } 

r[s, 

f . 

-> B°K*+] 

rW - 

* B* + K°] 

T[B+ - 

-»• B* + K*°] 

r[B+ - 

* B*°K+} 

r[s, 

+■ . 

-»• B*°K*+] 



v„ 

2 

(~) 2 ■ 



ri 


\aiJ 


This procedure can give the test for the factorization approach itself. 

The suppressed decays caused by the flavor changing neutral currents were studied in [29]. 

The CP-violation in the B c decays can be investigated in the same manner as made in the 
charged B decays. The expected CP-asymmetry of A(B± —> J/t^D ± ) is about 4 • 10 -3 , when 
the corresponding branching ratio is suppressed as 10 ~ 4 [30]. The reference-triangle ideology 
can by applied for the mo del-independent extraction of CKM-matrix angle 7 . However, the 
corresponding branchings are suppressed, e.g. Br(.£l7 —> D+D°) 10 5 . Thus, the direct 

study of CP-violation in the B c decays is practically difficult because of low relative yield of 
B c with respect to ordinary B mesons: <j(B c )/ct(B) ~ 10~ 3 . 

Another possibility is the lepton tagging of B s in the Bf —> B^^u decays for the study 
of mixing and CP-violation in the B s sector [31]. 
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Tabic 5: Branching ratios [16] of exclusive B+ decays at the fixed choice of factors: a\ = 1.20 
and a 2 = —0.317 in the non-leptonic decays of c quark, and a\ = 1.14 and a 2 = —0.20 in the 
non-leptonic decays of b quark. The lifetime of B c is appropriately normalized by t[B c ] ~ 0.45 
ps. 

kfode BR, % 

B+ -► J/i/jK + 0.011 

B c -> J/ifiK* + 0.022 

B+ -> D + W 0.0053 

B+ -> D+ZT 0 0.0075 

D+ -> ZT+D 0 0.0049 

D+ -> D*+D*° 0.033 

B+ -> D+D° 0.00048 

B+ -> DfTB° 0.00071 

B+ -> Z^+D 0 0.00045 

B+ -> D* s + T)*° 0.0026 

R+ -»• rj c Df 0.86 

5+ -> r/ c Ii: + 0.26 

5+ J/V^+ 0.17 

B+ J/^D; + 1.97 

B+ -»• r/ c Zt+ 0.032 

B+ -> 0.010 

R+ -> J/ , tpD + 0.009 

R+ -> 0.074 

R+ -> fi s °7r+ 16.4 

B+ -> fi s °p + 7.2 

B+ -> Rf 7T+ 6.5 

B+ -> R*°p+ 20.2 


We present here the current status of the R c meson decays. We have found that the various 
approaches: OPE, Potential models and QCD sum rules, result in the close estimates, while 
the SR as explored for the various heavy quark systems, lead to a smaller uncertainty due 
to quite an accurate knowledge of the heavy quark masses. So, summarizing, we expect that 
the dominant contribution to the B c lifetime is given by the charmed quark decays (~ 70%), 
while the 6-quark decays and the weak annihilation add about 20% and 10%, respectively. 
The coulomb-like aq/u-corrections play an essential role in the determination of exclusive form 
factors in the QCD SR. The form factors obey the relations dictated by the spin symmetry of 
NRQCD and HQET with quite a good accuracy expected. 

The predictions of QCD sum rules for the exclusive decays of B c are summarized in Table 
5 at the fixed values of factors «i 2 and lifetime. In addition to the decay channels with 
the heavy charmonium J/ip well detectable through its leptonic mode, one could expect a 
significant information on the dynamics of B c decays from the channels with single heavy 
mesons, if an experimental efficiency allows one to extract a signal from the cascade decays. 
An interesting opportunity is presented by the relations for the ratios in (10), which can shed 
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light to characteristics of the non-leptonic decays in the explicit form. 
We have found that the b decay to the doubly charmed states gives 


Br [B+ —> cccs] ~ 3.26%, 


so that in the absolute value of width it can be compared with the estimate of spectator decay 

[ 21 ], 


r[B 

r[s+ 


+ 


cc cs 


cccs 


SR 

spect. 


« 48 • 10 -15 GeV, 
« 90 • 10 -15 GeV, 


and we hnd the suppression factor of about 1/2. This result is in agreement with the estimate 
in OPE [21], where a strong dependence of negative term caused by the Pauli interference 
on the normalization scale of non-leptonic weak lagrangian was emphasized, so that at large 
scales one gets approximately the same suppression factor of 1/2, too. 

To the moment we certainly state that the accurate direct measurement of B c lifetime can 
provide us with the information on both the masses of charmed and beauty quarks and the 
normalization point of non-leptonic weak lagrangian in the B c decays (the cq and eq factors). 
The experimental study of semileptonic decays and the extraction of ratios for the form factors 
can test the spin symmetry derived in the NRQCD and HQET approaches and decrease the 
theoretical uncertainties in the corresponding theoretical evaluation of quark parameters, as 
well as the hadronic matrix elements, determined by the nonperturbative effects caused by the 
quark confinement. The measurement of branching fractions for the semileptonic and non- 
leptonic modes and their ratios can give information on the values of factorization parameters, 
which depend again on the normalization of non-leptonic weak lagrangian. The charmed quark 
counting in the B c decays is related to the overall contribution of b quark decays, as well as 
with the suppression of b —> ccs transition because of the destructive Pauli interference, which 
value depends on the nonperturbative parameters (roughly estimated, the leptonic constant) 
and non-leptonic weak lagrangian. 

Thus, the progress in measuring the B c lifetime and decays could improve the theoretical 
understanding of what really happens in the heavy quark decays at all. 

We point also to the papers, wherein some aspects of B c decays and spectroscopy were 
studied: non-leptonic decays in [32], polarization effects in the radiative leptonic decays [33], 
relativistic effects [34], spectroscopy in the systematic approach of potential nonrelativistic 
QCD in [35], nonperturbative effects in the semileptonic decays [36], exclusive and inclusive 
decays of B c states into the lepton pair and hadrons [37], rare decays in [38], the spectroscopy 
and radiative decays in [39]. 


6 B c production 

The (bc) system is a heavy quarkonium, i.e. it contains two heavy quarks. This determines 
the general features for the B c meson production in various interactions: 

1. Perturbative calculations for the hard associative production of two heavy pairs of cc 
and bb and 

2. A soft nonperturbative binding of nonrelativistic quarks in the color-singlet state can 
be described in the framework of potential models. 
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The two above conditions result in the suppression of the B c yield of the order of 10~ 3 
with respect to beauty hadrons. 

As was mentioned above, the consideration of mechanisms for the hadronic production 
of different spin P c -states is based on the factorization of hard parton production of heavy 
quarks ( bbcc) and soft coupling of (be) bound state [12]. In the first stage of description, the 
hard subprocess can be reliably calculated in the framework of QCD perturbation theory, 
while in the second stage the quark binding in the heavy quarkonium can be described in the 
nonrclativistic potential model assigned to the (6c)-pair rest system. The latter means that 
one performs the integration of the final quark state over the quarkonium wave function in 
the momentum space. Since the relative quark velocity inside the meson is close to zero, the 
perturbative matrix element can be expanded in series over the relative quark momentum, 
which is low in comparison with the quark masses determining the scale of virtualities and 
energies in the matrix element. In the leading approximation one considers only the first 
nonzero term of such expansion, so that for the 5-wave states the matrix element of the parton 
subprocess for the B c production is expressed through the perturbative matrix element for the 
production of four heavy quarks ( gg —> bbcc ) with the corresponding projection to the vector 
or pseudoscalar spin state of (6c)-system, which is the color singlet, and through the factor of 
radial wave function at the origin, R n s( 0), for the given quarkonium. The perturbative matrix 
element is calculated for the (bc) state, where the quarks move with the same velocity, i.e. one 
neglects the relative motion of 6 and c. 

For the P- wave states, the potential model gives the factor in the form of first derivative of 
the radial wave function at the origin, R' nL ( 0). In the perturbative part, one has to calculate 
the first derivative of the matrix element over the relative quark momentum at the point, 
where the velocities of quarks, entering the quarkonium, equal each to other. 

Thus, in addition to the heavy quark masses, the values of -R n s(0), R' nL ( 0) and a s are the 
parameters of calculation for the partonic production of B c meson. In calculations we use the 
wave function parameters equal to the values shown in table X and R' 2 p(0) = 0.50 GeV 3 / 2 . 
The value of R( 0) can be related with the lcptonic constant, /, so that we have 

fis = 0.47 GeV, f 2S = 0.32 GeV. and f n = J-^-R nS ( 0) . 

V nM n 

In the approximation of the weak quark binding inside the meson one has Mp c = mb + m c , 
so that the performablc phase space in calculations is close to physical one at the choices of 
mb = 4.8 GeV, m c = 1.5 GeV for the 15-state, nib = 5.1 GeV, m c = 1.8 GeV for the 25-state, 
mb = 5.0 GeV, m c = 1.7 GeV for the 2P-state. 

At large transverse momenta of the B c meson, pp 3> Mb c , the production mechanism 
enters the regime of 6-quark fragmentation (see fig. 1), so that the scale determining the QCD 
coupling constant in hard 66 production is given by /i|. ~ M\ + Pp, and in the hard fragmen¬ 
tation production of the additional pair of heavy quarks cc we get p 5c ~ ?7V This scale choice 
is caused by the high order corrections of perturbation theory to the hard gluon propagators, 
where the summing of logarithms over the virtualities leads to the pointed // values. Therefore, 
the normalization of matrix element is determined by the value of ay(/q; fe )ay(/x gc ) ~ 0.18 • 0.28. 
In calculations we use the single combined value of a s = 0.22. 

The parton subprocess of gluon-gluon fusion gg —> B+ + 6 + c dominates in the hadron- 
hadron production of B c mesons. In the leading approximation of QCD perturbation theory 
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it requires the calculation of 36 diagrams in the fourth order over the a s coupling constant. 

By the general theorem on factorization, it is clear that at high transverse momenta the 
fragmentation of the heavier quark Q —> ( Qq ) + q must dominate. It is described by the 
factorized formula 


da 

dp T 


dd{m gg 


dki 


— \kr-rr,* ■ DV-IWfx;»)-, 

oc 


( 11 ) 


where p is the factorization scale, da/dkx is the cross-section for the gluon-gluon production 
of quarks Q + Q, D is the fragmentation function. 

The calculation for the complete set of diagrams of the 0(a£)-contribution [12] allows one 
to determine a value of the transverse momentum p™ m , which is the low boundary of the region 
where the subprocess of gluon-gluon Z? c -meson production enters the regime of factorization for 
the hard production of 66-pair and the subsequent fragmentation of 6-quark into the bound 
(6c)-state, as it follows from the theorem on the factorization of the hard processes in the 
perturbative QCD. 

The p r rf nn value turns out to be much greater than the Mb c mass, so that the dominant con¬ 
tribution into the total cross-section of gluon-gluon H c -pro duct ion is given by the diagrams of 
nonfragmentational type, i.e. by the recombination of heavy quarks. Furthermore, the convo¬ 
lution of the parton cross-section with the gluon distributions inside the initial hadrons leads to 
the suppression of contributions at large transverse momenta, as well as the subprocesses with 
large energy in the system of parton mass centre, so that the main contribution into the total 
cross-section of hadronic _B c -pro duct ion is given by the region of energies less or comparable 
to the -B c -meson mass, where the fragmentation model can not be applied by its construction. 
Therefore, one must perform the calculations with the account for all contributions in the 
given order under consideration in the region near the threshold. 

The large numeric value of p“ in means that the majority of events of the hadronic B^- 
production does not certainly allow the description in the framework of the fragmentation 
model. This conclusion looks more evident, if one considers the _B c -meson spectrum over the 
energy. 

The basic part of events for the gluon-gluon production of B c is accumulated in the region 
of low z close to 0, where the recombination dominates. One can draw the conclusion on 
the essential destructive interference in the region of z close to 1 and px < p™ 111 , for the 
pseudoscalar state. 

We have considered in detail the contributions of each diagram in the region of 3 —> 1. In 
the covariant Feynman gauge the diagrams of the gluon-gluon production of Q + Q with the 
subsequent Q —> (Qq) fragmentation dominate, as well as the diagrams when the qq pair is 
produced in the region of the initial gluon splitting. However, the contribution of the latter 
diagrams leads to the destructive interference with the fragmentation amplitude, and this 
results in the ’’reduction” of the production cross-section in the region of z close to 1. In the 
axial gauge with the vector nh = pQ this effect of the interference still manifests itself brighter, 
since the diagrams like the splitting of gluons dominate by several orders of magnitude over the 
fragmentation, but the destructive interference results in the cancellation of such extremely 
large contributions. This interference is caused by the nonabelian nature of QCD, i.e., by the 
presence of the gluon self-action vertices. 

The using of CTEQ5L parameterization for the structure functions of nucleon [40] leads to 
the total hadronic cross-sections for the B c mesons of about 0.8 pb that accepts contributions 
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-3 



GeV 

Figure 3: The differential cross-section for the meson production in gluon-gluon collisions 
as calculated in the perturbative QCD over the complete set of diagrams in the 0(a*) order 
at 200 GeV. The dashed and solid histograms present the pseudoscalar and vector states, 
respectively, in comparison with the corresponding results of fragmentation model shown by 
the smooth curves. 


15 






from: 


1S Q IS! 2 S 0 2 Si 

0.19 fib 0.47 /ib 0.05 /ib 0.11 /ib 

After the summing over the different spin states, the total cross-sections for the production 
of P-wave levels is equal to 7% of the 5-state cross-section. 

At LHC with the luminosity C = 10 34 cm _ 2 s _1 and y/s = 14 TeV one could expect 
4.5 • 10 10 P+ events per year. 

Nevertheless, the P-wave states could be of a particular interest due to their radiative 
decays with relatively energetic photons (around 500 MeV in the B c rest frame). For P-wave 
states, the leading color-singlet matrix element and the leading color-octet matrix elements 
are both suppressed by a factor of v 2 (relative velocity of the charmed quark) relative to the 
color-singlet matrix element for S-wave that can enhance the P-vawe contribution. 

In the Figure 4 the da/dy distribution is presented. The y distribution shows the maximum 
in the central region however, considereing the experimental observability of the B c states, one 
should care about the momentum of the meson to ensure the reasonable 7 -factor for B c and 
visible separation of B c decay vertex from the primary one. 

The right part of the same Figure shows the dependence of the B c momentum on the angle. 
One can see that the central region is dominated by low-momentum mesons, which makes the 
observation of these states in the central region quite difficult task. On the other hand, the 
forward-backward (LHCb acceptance) regions are dominated by very energetic mesons. 

In the region of the LHCb acceptance (O < 17°) the expected number of events with 
the nominal luminosity of 2 x 10 32 is about 10 9 per year. Taking the value of ~ 0.1 as an 
approximation for the reconstruction efficiency of the decay B c —> J/T tt —> /i/i 7 r, for example, 
one gets the total amount of 2 x 10 4 reconstructed events per year. 

The topology of the events with P c -mesons production is somewhat specific due to extreme 
kinematics which, particularly, is responsible for the enhancement of the forward-backward 
regions. The main feature of these events topology consists in the strong corellation in the 
direction of B c and associated D- and 6 -mesons momenta. One should expect an associated 
production of all three heavy mesons in the same hemisphere and, moreover, in sufficiently 
narrow cone. 


7 Conclusions 

The family of (be) mesons contains 16 narrow states. The S-wave ones will be produced in 
pp collisions at LHC energies with relatively large cross-sections, ~ 0.1 /ib. The total cross- 
section of the B c production, taking into accout the cascade decays of the narrow excited 
states, can be as high as ~ 1 /ib. This value is more than order of magnitude larger than at 
Tevatron energy. 

With the total luminosity of about C = 10 34 cm _ 2 sec _1 , one could expect the total amount 
of B c mesons produced of the order of Nb c ~ 5 • 10 10 per year. 

The forward-backward regions of the P c -mesons production are more favourable in the 
view of experimental observation due to the strong Lorentz boost of the initial parton system. 
One can expect ~ 10 9 B c events per year inside LHCb acceptance. This amount is quite 
sufficient to study the spectroscopy and various decay modes, as well as the life-time of the 
ground state. 
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Figure 4: da/dy (left) and B c momentum vs. angle (right) for the _B+(l£o) at LffC energy. 

The inclusive decay mode B c —> J/TX has a branching of about 17%, in comparison with 
~1% in Bd s -mesons decays. These channels produce a very well visible signatures and one 
could expect ~ (2 — 4) • 10 4 events in the decay mode B c —> J/T7T —*• /lyi r and ~ 8 • 10 4 events 
in the B c — > J/Tpz/ —> jip/iv mode. The Cabbibo-suppressed mode B c — > .J/'i’K —> pyK can 
be observed at the level of ~ 10 3 events. 

The most probable 6-spectator decay modes are saturated by the two-body decays B c —> 
B^^ and B c —> B^ p ± , that makes these channels quite interesting for the study of the 
c-quark decays in the £> c -meson. The estimated yield of the reconstructed events could be 
~ 10 3 events per year in the B c —> B s 7r ± mode. Approximately the same amount of events 
can be observed in other modes with B s in the final state, although with worse background 
conditions. 

Rare decay modes, e.g. B c —» D ± D° , could be interesting in view of CP-violation studies, 
however the preliminary estimations of the detection efficiencies and branchings involved are 
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not too optimistic and additional studies of the reconstruction and selection procedures are 
required. 

Authors are very grateful fot the contribution of Alexandr Berezhnoj to this work. We 
also would like to thank LHCb collaboration for warm hospitality and IHEP group (Protvino) 
in LHCb that provides a support for this work. 
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